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POWER-SERIES  SOLUTIONS  FOR  FLOWS  OF  AN 


IDEAL  DISSOCIATING  GAS 


1 .  Introduction 

This  report  deals  with  the  solution  of  certain  problems 
in  fluid  mechanics  by  power-series  expansion  of  the  solution 
in  the  Independent  variable (s) .  The  method  is  directly 
related  to  the  well-known  Frobenius  method  for  determining 
analytic  solutions  of  linear  ordinary  differential  equations 
(see,  e.g.,  Agnew  I960).  However,  hero  we  apply  it  to  nonlinear* 
systems  of  differential  equations  in  as  many  as  three  Independent 
variables.  In  consequence ,  the  recursion  formulas  for  the 
series  coefficients  are  relatively  complicated,  and  an 
electronic  computer  is  required  to  affect  and  to  store  their 
solution .  . 

There  are  some  texts  on  numerical  methods  which  discuss 
power-series  methods  for  the  solution  of  nonlinear  ordinary 
differential  equations  (e.g*,  Henrici  1964),  but  usually 
only  as  .a  prelude  to  the  description  of  more  "practical 
methods.  Most  applications  of  series  expansions  to  nonlinear 
problems  have  been  fox  the  purpose  of  clarifying  some  Local 
behavior  of  the  solutions  and,  thus  requiring  only  the  first 
few  terms  of  the  series,  have  been  carried  out  analytically 
rather  than  numerically.  Nevertheless,  the  use  of  series 
methods  in  fluid  mechanics  has  a  considerable  history. 

Most  of  this  history  is  related  to  the  problem  of  the 
steady  supersonic  flow  of  a  perfect  gas  past  a  blunt  body. 
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On  prescribing  the  bow  shock  shape  (analytically) ,  one  can 
substitute  approximate  aeries  expansions  for  the  flow  variables 
into  the  governing  equations  and  determine  their  coefficients 
recursive!'/ .  .  nopea  v.aa  rubies)  in  which  the  series 

converges  encompasses  the  body  contour  (which  must  b 
determined  as  part  of  the  solution)  and  the  sonic  line,  beyond 
which  this  ic  .uc:.nr.'  cam  be  continued  by  the  method  of  character¬ 
istics.  First  to  suggest  ciue  approach  ware  Lin  and  Rubincv 
(194 H).,  Early  attempts  uo  implement  it  were  made  by  Dugutvdji 
(1948) ,  Lin  and  Shell  (1951),  and  Cabannec  (IS f  .495C-  ? 

The  first  to  use  a  computer  hr:  able  to  generate  a 

substantial  number  or  series  coefficients  was  Richtmyer  (1357) . 

,*  7;,:s  Dyke  (1953)  soon  gave  evidence  tL-r.r.  fcfco 
series  do  not  converge  Lccy.  Bh  ascribed  this 

failure  to  the  appearance  in  the  analytic  continuation  >r?. 
the  series  upstream  of  the  shock  a  singularity  which  is  closer 
to  the  center  ox  che  expansion  than  is  ths  ,  subsequent 
studi.«ss  were  therefore  directed  at  extending  the  rag;  :  of 
convergence  of  the  series  so  that  they  would  have 
pr  V--:.'  h"!./..  .  .\>SV ,•  I9T5.)  .  Van  Dyke  (l?6fc) , 

u<.:, ■.  U365)  ,  Leavitt  (1963),  and  Moran  (1970)  all  devised 
successful  eolu lions  of  this  probl: 

The  .  •  noted  above  are  all  based  on  expansion  oi 

me  flow  properties  m  series  in  each  of  the  independent 
variables.  In  1908,  Bint  its  suggested  that  the  growth  of  ;■>. 
two-dimensional  laminar  boundary  layer  an  a  pressure  gradient 
could  be  determined  by  expanding  cue  solution  in  a  pow, 
series  rn  the  variable  runnii.g  along  the  body.  The  coefficient 


would  depend  on  the  coordinate  normal  to  the  body  and  would  be 
determined  by  numerical  solution  of  appropriate  ordinary 
differential  equations  (see  Schlicting  1958  for  subsequent 
developments  of  this  '’Blasius  series”  approach)  ,  Van  Dyke 
(1965)  applied  a  similar  partial-expansion  method  to  the 
blunt-body  problem  with  excellent  results.  He  also  used  this 
procedure,  which  he  calls  the  series  truncation  method,  on  the 
Navier-Stokes  equations  for  the  flow  past  a  circular  cylinder 
(Van  Dyko  1964,  Underwood  1969) .  Davis  (1967)  studied  flow 
past  a  flat  plate  with  the  same  method.  In  the  meantime, 

Conti  and  Van  Dyke  (1966,  1969)  have  treated  inviscid 
flows  of  dissociating  gasses  past  blunt  bodies  with  partial- 
expansion  methods,  while  Moran  and  Van  Moorhem  (1969,  1970) 
used  a  full-expansion  method  for  the  three-dimensional 
(two  space,  one  time)  flows  of  a  perfect  inviscid  gas  which 
are  generated  when  a  plane  shock  wave  strikes  either  a 
stationary  blunt  body  or  one  that  is  in  steady  supersonic 
flow. 

The  above  listed  flow  situations  have  in  common  a  smooth 
variation  of  the  flow  properties  with  the  variables  in  which 
they  are  expanded.  To  be  sure,  boundary  layers,  in  which  the 
properties  vary  rapidly,  do  develop  in  the  flow  of  an  inviscid 
dissociating  gas  paBt  a  blunt  body  and  in  the  high-Reynolds 
number  flow  of  a  viscous  fluid  past  a  solid  body.  However, 
the  expansions  used  by  Conti  and  Van  Dyko  (i.966,  1969), 

Davis  (1967),  and  Underwood  US69)  in  treating  these  problems 
were  only  in  the  variables  running  parallel  to  these  layers. 

On  the  other  hand,  while  Conti,  Van  Dyke  and  Davis  were  able 
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to  describe  their  boundary  layers  quite  successfully,  the 
series  of  Underwood  converge  only  in  the  low-Reynolda -number 
regime,  where  the  boundary  layers  are  relatively  thick  and 
the  property  variations  relatively  smooth. 

The  present  investigations  were  directed  at  elucidating 
the  ability  of  the  power-series  methods  to  treat  problems  with 
embedded  boundary  layers  (as  distinguished  from  cases  in  which 
coordinate  stretching  has  narrowed  the  focus  to  the  boundary 
layer  itself,  as  in  Blasius's  (1908)  problem).  The  problems 
selected  involved  cue  iibri.ua;  flow  of  a  dissociating 

gas  (specifically,  the  "ideal  dissociating  gas"  of  Lighthill 
(19S7)  and  Freeman  (1958)).  Behind  a  strong  shock,  the  finite- 
rate  relaxation  of  such  a  gas  towards  equilibrium  creates  a 
region  in  which  the  f lo *  properties  may  vary  rapidly  relative 
to  their  variations  elsewhere  in  the  flow.  That  is,  trie 
characteristic  relaxation  length  may  be  small  compared  to  the 
characteristic  geometric  length.  Boundary  layers  also  appear 
near  stagnation  points  when  the  ratio  cf  these  characteristic 
lengths  is  at  the  other  extreme. 

The  equations  which  govern  all  the  flows  studied  ere  set 
down  in  Section  2.  To  illustrate  the  basic  techniques  we 
employ,  the  relatively  simple  problem  of  determining  tha  flow 
properties  behind  a  strong  plane  shock  is  examined  in  detail 
m  tlis  next  three  Sections .  First  the  recursion  forim.’.l.. ; 
for  the  series  coefficients  are  derived,  then  the  nanipulaticns 
required  to  mike  maximum  uae  of  the  coefficients  determined  are 
described,  \nd  finally  the  results  obtained  are  examined. 


4 


Sections  6  and  7  are  devoted  to  two-dimensional  steady  tlows 
past  pointed  and  blunt  bodies.  Finally,  in  Section  8,  we 
treat  the  three-dimensional  (two  space,  one  time)  flow  which 
results  when  a  plane  shock  reflects  from  a  blunt  body. 

It  would  appear  from  our  results  that  series-expansion 
methods  are.  of  somewhat  limited  utility  in  describing  flows 
with  imbedded  boundary  layers.  While  a  full-expansion 
approach  is  generally  easier  to  program,  it  may  be  necessary 
to  use  a  partial  expansion,  with  no  expansion  in  the  variable 
running  normal  to  the  layer.  This  is  certainly  so  if 
singularities  appear  in  the  flow  field,  and  even  then  what 
amounts  to  a  local  analysis  of  the  singularity  (cf.  Conti 
and  Van  Dyx.e  1966 1  1969)  is  necessary  to  make  the  method 
useful.  On  the  other  hand,  where  the  series  can  be  made  to 
converge,  their  accuracy  is  superb,  this  justifying  much  of 
the  effort,  required. 
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2 .  Basic  Equations 


The  forms  of  the  conservation  equations  convenient  for  the 
present  investigations  are  as  follows  (Vincenti  and  Kruger 
1965,  p.  246): 


g|r  +  p  div  v  *  0  (1) 

Dv 

PDt  f  grad  P  “  0  (2) 

-H  •  81  - 0 

Here  p,  p,  h,  v  are  the  fluid  pressure,  density,  enthalpy, 
and  velocity,  and  D/Dt  is  the  convective  derivative. 

b't  "  H  +  grad  (4) 

For  Lighthill *s  ideal  dissociating  gas,  the  species  conservation 
equation  is 

25.  *  5L*Z£L  (5) 

Dt  T 

* 

where  a  is  the  degree  of  dissociation ,  a  its  equilibrium 
value,  and  t  a  relaxation  time.  The  law  of  mass  action  becoros, 
Lighthill  '  a  (1957)  approxi.mat ?  ^n, 


in 
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a*2 


ei  ,-vt 


(6) 


where  T  is  the  temperature  and  ,  0^  are  constants  of  the 
gas,  Freeman's  (1958)  equation  for  the  relaxation  time  .nay 
be  written 


t  «* 


l-o* 


Cp5T*  l-(l-a*){l-a) 


(7) 


where  C  Is  e  constant.  Finally,  the  equations  of  state  are, 
with.  R  the  gas  constant, 


p  “  p  (l+o }  RT 

(8) 

h  -  (4+o)  RT  +  o8d 

(9) 

the  latter  of  which  also  being  an  approximation  due  to  Lighthill 

(19571  . 

Across  shock  waves,  the  conservatiort  equations  become  the 

jump  conditions 

IP  (vR  -  x4U  *  0 

(1*-) 

(P  ♦  p(vft  -  x#)2]  •  0 

Ul> 

tv.]  -  0 

(12) 

7 


3 .  Relaxation  Behind  a  Steady  Normal  Shocks  Development  o' 
Series  Solution 

The  simplest  case  we  have  considered  is  the  one -dimens iontl 
flow  thru  a  standing  normal  shock  wave.  Let  x  be  the  distance 
downstream  of  the  shock.  With  D/Dt  «  vd/dxf  the  conservation 
equations  (l)-(3),  together  with  the  jump  conditions  (10),  (11), 
and  (13) ,  may  be  integrated  to  yield 


Pv  “ 


P  +  P«oVa.V  +  P.  +  PWVco‘ 
h  +  v2/2  »  hw  +  v„2/2 


where  the  subscript  *  denotes  values  upstream  of  the  shock. 
The  species  conservation  equation  (5) ,  however,  remains  in 
differential  form; 


da  _  a* -a 


*  CT  [  (p-pa)e  d 


-e./T  _  p2a2 


in  which  we  have  incorporated  equations  (6;  and  (7) .  From  the 
jump  condition  (14) ,  we  get  the  initial  condition 


a  (0)  *  a 


The  series  solution  of  these  equations  begins  with  the 
assumption  that  each  of  the  variables  may  be  expanded  in  a 
Taylor  series  of  the  form 


A  (x) 


iL 4  A.X^ 

j"0  p 
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We  then  obtain  recursion  formulas  for  the  coefficients  by 

substituting  the  series  into  the  governing  equations  and 

v 

collecting  terms  of  like  order  in  x  .  Prom  the  linear  equation 
(16),  this  yields 

Pm  +  p«v>»vm  m  <P-  +  p-v-2)  «„0  (21) 

with  no  trouble  whatsoever.  Equations  with  terms  quadratic  in 
the  unknowns  are  only  a  little  more  difficult.  Prom  (9),  (15), 
and  (17)  we  obtain 

hm  ’  <ETm  +  °iVi  +  9d°m  (22) 


m 

Lo .  v  .  »  o  v  6  - 
q  m-i  «*>  mO 


(23) 


ra 


£ 

i=0 


lv.v  < 
2*  i  m-i 


*  +  1  v  )  <5 


mO 


(24) 


However,  equations  like  (8),  with  &.  triple  product  of  unknowns, 
and  the  highly  nonlinear  (18)  ,  are  rather  nasty  to  work  with. 

We  generally  find  it  worthwhile  to  introduce  enough  auxiliary 
variables  so  that  no  equation  contains  terms  more  nonlinear 
than  quadratic  before  expanding  the  variables  therein  in 
Taylor  series.  Thus,  with 

b  5  pa  (25) 


eq  ua 1 1 on  ( 3 )  become s 
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P  *  (P  +  n)RT 


(26) 


and  thoir  Taylor-serias  expansions  yield  the  recursion  formulas 


m 

bn>  *  Z  piVl  (27) 

i«0 

m 

p»  -  *  £  ‘pi +  bi>Vi  <2e) 

i*0 

To  simplify  (18) ,  we  further  define 

f  £  6~VT  (29) 

g  5  0d/T  (30) 

w  =  T6  (31) 

d  2  (p-b) f-b2/pd  (32) 

Then  (18)  becomes 

-  Cwd  (33) 

Equations  (30),  (32),  and  (33)  are  easily  expanded  in  series  to 
give 

m 

£  «iVi  •  6Ao  (34) 

m 

-  biVi/p<?  <35> 

i-0 

m  m 

ia.v  ,  *  C  L  w.  -<i  i  (36) 

i«0  1  ra~i  i-1  i~1^n"i 
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To  help  expand  (29)  and  (31),  we  follow  Leavitt  (1366)  and  first 
differentiate  with  respect  to  x: 


df 


(37) 


dw  _  aw  dT 
Hx  6T  53T 


(38) 


These  differential  equations,  on  series  expansion,  yield  the 
recursion  formulas 
m 

»*»  *  -  E  lgi  Vi  <»> 

i=0 


m 


m 


i«l 


xw. 

l 


m-i 


6  L 

i=»l 


iT<Wm-4 

i  m-x 


(40) 


Note  that  equations  (36)  ,  (39)  ,  and  (40)  are  valid,  only  for 
m  >  0.  For  m  °  0,  they  are  replaced  by  the  initial  conditions 


an  =  a 
0  00 


fQ  -  e~‘3d/T0 


w0  "  V 


(36a) 

(39a) 


(40a) 


which  are  derived  from  "  qua  t  A  or  (19)  and  (after  setting  xmQ) 
from  the  definitions  (29)  and  (31) . 

A  pr  oblem  often  more  dii'ficuJ  t  Loan  the  formation  of  the 
recursion  formulas  amonst  the  series  coefficients  is  the 
do  termination  of  rhe  order  in  which  they  may  be  solved.  A 


•- 1  tuv  a.  V«.<, 

first 

p  tovrar 

do  :  t. 

Oi  nxyr; 

esfc  0 

rdar 

in 

h  cl: 

on#  hop 

hh  (a 

•s  in 

tfc 

T>  £>'  *» 

oh 

ao a file 

iencs 

X  j’ '  it 

asqtfisno 

e  acc 

reveal a 

what 

oqt 

at! 

on  3 

iftC&t 

Table  .1 

,  wa 

aoa 

rha 

t  eg 

u&ti  0 

o  f  car  v/hich 

ptr/ 

V 

V 

m 

simultaneous  solution  cf  e^ue.tiona  (31)  f  (22)  ,  £23)  »  (24)  / 

(27)  ,  and  (25) .  Tfoesa  5,.  c&«  fea  found  from  (34) ,  fw  from  (39) 
or  (39a) ,  <3  from  (35) ,  and  from  (40)  or  (40a).  This 
routine  is  fallowed,  for  m,  »  0,  1*  7,  . .  * »  end  thus  yield* 
coefficients  of  the  neriso  expans  ions  of  arbitrarily  hieh  orde.c. 

The  equations  to  be  solved  too-  all  linear/  except  wheu 
m  *•  0,  and  even  then  an  analytic  solution.  is  not  hard  to  obtain. 
However,  the  coefficients  of  the  linear  equations  rapidly 
increase  in  complexity  with  incres.sjlarcf  :a.  r?hu.3r  if  ore  wishes  to 
obtain  more  than  the  first  £tv.v  reri.r:-*  coefficients,  he  raust 
use  a  digital  computer .  '°x ograuircrq  the  rj*cur.->ion,  rt iutiens  rs 

straightforward „  Ls:-.vi to  £X3-t'$>  ana  <f'i.a3y.v  ••‘.s.rcilvy  the 
frocjrarfl»iin<j  by  writing  snlroufintr  hr;  r>vtc5a:-v.rn<y  tha  tert  of 
order  :c  from  series  wttct  are  >•.{••::  dnnbio  c.r  t;  i;;  .w::.  product  or 
exponential,  etc...  o.  oth-u  snri***.  howewy  do  n-r.i  find  the 
uiroet  y  rogr  nr  a:  ti  -.ho  r«.cn*..r£ic^  fcrmif  as  to  entrees  on  to  ba 

worth  the  pr->  »n  e icreg--*  and  cot nrv.no *rr>r.  1 '  ,-c-  ‘  r.c  rr. rf,  pry  for 


Utilization  of  Sorias  Solution 
Once  the  series  coefficients  are  known,  the  next  step  is 
to  use  the  series  to  calculate  the  flow  properties  at  points 
of  interest.  Where  the  series  converge,  they  may  be  used 
directly,  and  the  accuracy  of  the  results  may  be  estimated  by 
examining  the  effect  of  varying  the  number  of  terms  after  which 
the  series  is  truncated. 

Unfortunately,  the  series  we  encounter  typically  have 
disappointingly  small  regions  of  convergence .  The  singularities 
which  limit  convergence  are  usually  outside  the  domain  in  which 
the  series  is  of  physical  interest.*  The  one  which  is  troublesome 
in  the  relaxation  problem  appears  (from  examination  of  the 
series),  to  be  located  at  some  negative  value  of  x,  which  would 
put  it  upstream  of  the  shock,  where  the  flow  is  uniform  and  not 
described  by  the  series  at  ail.  This  artificiality  does  not 
make  the  singularity  any  less  dangerous  ;■  the  region  in  which 
a  series  converges  depends  only  on  the  distance  from  the  center 
of  the  expansion  to  the  nearest  singularity,  and  is  otherwise 
independent  of  the  location  of  that  singularity. 

Thus  r  for  the  series  solution  tc  be  useful,  vze  must  be  able 
to  continue  it  analytically  beyond  its  ba3ic  region  of  conver¬ 
gence,  given  (numerically)  only  the  first  several  coefficients  of 
the  series.  Lewis  (1965)  studied  a  number  of  approaches  to  this 


problem.  Probably  the.  most  powerful  is  simply  to  introduce  a 


new  independent  variable 


so  that,  when  the  series  is  reexpanded 


*When  they  are  within  that  domain,  our  method  generally  fails; 
see  Sections  6  'and  7  below. 
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in  the  new  variable,  there  is  a  more  favorable  relationship 
between  the  locations  of  the  singularities  of  the  function  and 
the  domain  of  interest.  Such  transformations  were  developed 
further  by  Leavitt  (1968)  ,  who  made,  use  of  Domb's  (1965) 
procedure  for  locating  and  classifying  the  singularity  closest 
the  origin  of  a  function  known  only  by  the  coefficients  of  its 
power-series  expansion. 

Lewis  and  Leavitt,  both  obtain  splendid  results  for  the 
supersonic  flow  of  a  perfect  gas  past  a  blunt  body.  The  pro¬ 
cedures  we  have  employed,  however,  are  based  on  the  theory  of 
continued  fractions.  Given  the  coefficients  ot  a  power  series 

2 

ao  +  alx  +  a2x  +  •••  53  (41) 

we  can,  under  certain  restrictions  cn  the  an's,  find  numbers 
an  such  that  the  power-series  expansion  of  the  continued 
fraction 

a0/(l  +  a^x/d  +  a,x/  (1  +  a3x/(l  +  . (42) 

agrees  with  the  left  side  of  (41) .  hn  efficient  procedure  for 
so  doing  is  the  "quotient -difference  algorithm"  (Henrici  1963)  , 

The  theory  of  continued  fractions  is  full  of  half -promises 
concerning  the  efficacy  o:  representing  functions  by  expansions 
like  (42).  It  is  known,  for  example,  what  conditions  on  the 
coefficients  of  the  pever-seriee  expansion  of  f(x)  are 
necessary  c  id  sufficient,  t ->r  it  to  possess  a  continued-fraction 
expansion  like  (42)  (Wall  1948,  Chap.  XI).  However,  these 
conditions  take  the  form  that  certain  matrices  whose  elements 
are  the  power-series  coefficients  must  be  noneingular.  Such 


15 


conditions  cannot  be  verified  when,  as  in  our  case,  the 
coefficients  are  known  only  numerically,  because  of  roundoff 
errors.  It  is  also  known  (Wall  1948,  p.  399)  that  every  function 
analytic  at  the  origin  nas  a  continued-fraction  expansion  oZ 
a  form  somewhat  more  general  than  (42);  specifically,  the  quantities 
a^x  must  be  replaced  by  a^x®k,  where.  8^  is  a  positive  integer  not 
necessarily  equal  to  unity.  Eowever,  the  numerical  difficulties 

required  to  determine  the  '3  also  appear  to  be  insurmountable. 

k 

Thus,  all  we  can  do  is  to  assume  that  the  power  series  with  which 
we  deal  have  continued-fraction  expansions  of  the  form  (42) . 

The  convergence  properties  of  such  expansions  are  rather 
spectacular  in  certain  cases,  /hi ch  of  course  is  what  motivates 
us  tu  use  the  expansions.  If  a  function  f (x)  ran  ba  expanded 
in  the  form  (42)  ,  and  if  the  coefficients  a^  approach  a  finite 
limit  ,  a  (sav)  it  cCTn  be  show  (Wall  1948,  Chap.  XI)  that  the 

<XJ  ' 

expansion  converges  to  f  (x)  eve rywhere  in  the  complex  x-plane 
except  at  the  poles  of  the  f unrf.ons  (if  any)  and  along  a  branch 
cut  which  begins  at  the  point  x  =  -l/4a^  and  proceeds  to  infinity 
along  a  ray  from  the  origin.  Ho.  e-er,  it  is  also  possible  to 
constrv  .t  continued- fraction  nxpannions  which  diverge  within 
the  circle  of  convergence  of  the  corresponding  power-series 
expansion.  Fortunately,  such  examlcs  appear  to  be  pathological. 

We  have  yet  to  encounter  i  power  s-;  ri-.s  which  converged  met  ter 
than  tne  cor;  esponding  cont  imwad  /  ract  \or. . 


5.  Relaxation  Dehind  a  Steady  Normal  Shock: _ Results 

Calculations  have  been  carried  out  as  described  above  for 
the  three  cases  identified  in  Table  2,  which  were  previously 
studied  by  Freeman  (1958)  with  a  more  conventional  numerical 
method.  Before  carrying  out  the  computations,  the  equations  were 
of  course  made  dimensionless.  In  particular,  the  independent 
variable  in  which  the  series  are  amended  is  not  the  physical 
distance  downstream  of  the  shock  x,  but 


The  first  twenty-five  coefficients  of  the  power-series 
expansion  of  the  temperature  behind  the  shock  for  the  first 
case  of  Table  2  are  listed  in  Table  3,  along  with  the  coefficients 
of  the  corresponding  continued-fraction  expansion.  Note  that 
the  latter  tend  to  remain  of  the  same  order  of  magnitude, 
which  is  something  of  a  convenience.  However  (unless  the  approach 
is  completely  obscured  by  roundoff  errors) ,they  evidently  do 
not  approach  a  limit,  so  that  tne  theorems  quoted  in  the  previous 
section  shed  no  light  on  the  convergence  properties  of  the  con¬ 
tinued-fraction  expansion. 

Of  course,  empirical  evidence  cn  the  utility  of  both  the 
power-series  and  continued-fraction  expansions  may  be  obtained  by 
varying  the  number  of  terms  used.  Data  relevant  to  the 
expansions  whooe  coefficients  are  tabulated  in  Table  3  are 
contained  in  Tables  and  5.  It  is  clear  that  the  power  series 
is  useless  even  for  x*  »  8.0  ,  although  examination  of  the 
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coefficients  by  Domb'p  (1965)  procedures  indicates  that  the 
series  should  converge  up  to  about  x*  *  11.5.  The  continued- 
fraction  expansion,  on  the  other  hand,  is  quite  satisfactory  for 
x*  well  above  100.  Of  course,  !*■«  accuracy  eventually  diminishes 
with  increasing  x*,  due  tc  roundoff  errors  if  nothing  else; 
the  results  listed  certainly  appear  to  be  trying  to  converge. 

On  the  other  hand,  the  fact  that  the  continued-fraction  results 
appear  to  oscillate  with  diminishing  amplitude  aB  the  number 
of  terms  is  increased  is  somewhat  misleading.  As  can  be  seen 
.rom  Table  5,  the  exact  result,  (obtained  by  what  is  essentially 
Freeman's  (1959)  method  cf  direct  integration)  does  not 
necessarily  fall  within  the  range  of  the  oscillations.  This, 
again,  is  probably  due  to  roundoff  errors  in  the  continued- 
fraction  coefficients  and/or  in  the  evaluation  of  the  expansion, 
■Results  obtained  from  the  continued-fraction  expansion  for 
T  are  compared  with  the  numerically  exact  results  for  cases  1 
and  2  in  figure  1  (the  results  for  cases  2  and  3  are  much  the 
same) .  To  give  some  indication  of  the  apparent  convergence 
of  the  expansion,  results  obtained  both  by  using  34  terms  and 
by  using  35  terms  are  plotted,  the  difference  between  them  being 
indicated  by  a  bar.  Tc  is  seen  that  the  results  benave  as 
described  above  in  both  of  cases  .1  and  2,  except  that  the  range 
of  x*  for  which  the  continued  -fraction  expansion  is  useful 
is  more  restricted  i...  .;.uc  ?  [dad  m  case  3}  than  in  case  1  . 
However,  the  power -serxes  expansions  in  the  last  two  cases  do 
note  converge  beyond  x'*  1.5,  whicn  is  far  below  the  largest, 

x*  for  which  tiie  cot  re  noon  sing  continued  fraction  io  useful, 
result;-,  for  the  other  properties  behind  the  shock  are  quite 
sin  i.  far  . 

IS 


Of  course,  the  series-expansion  method  is  no;  really  practical 
for  the  problem  of  finding  the  flow  behind  a  shock,  in  a  relaxing 
gas;  it  cannot  compete  with  ordinary  numerical  integration 
either  in  range  of  utility  or  in  time  (the  present  method  took 
about  4  seconds  on  a  CDC  6600  to  obtain  results  equivelent  to 
those  obtained  hy  a  Runge-Kutta  method  in  2  seconds.}  However, 
the  series  method  is  much  more  attractive  in  two-  and  t  .ree- 
dimansional  problems,  for  which  alternative  procedures  are  at 
least  as  costly  and  much  less  accurate.  For  the  presen;,  the 
main  conclusion  we  can  draw  is  that  the  series  method  can  be 
useful  through  the  most  interesting  part  of  the  shock  layer. 

The  failure  of  the  method  to  describe  the  solution  accurately 
as  it  approaches  equilibrium  may  simply  be  due  to  the  larg» 
arguments  involved  in  the  expansion. 


6 .  Steady  Supersonic  Flow  Past  vtedges  arc  Cones 

The  next  "lost  difficult  problsss  we  attempted  was  the  steady 
flow  behind  a  shock  attached  to  a  wedge  or  cone.  In  spherical 
coordinates,  equations  (l)-{5)  reduce  to 


3(pu)  .  {1  +e)  ou  .  pv  , 1  «Pv 
lx~  +  “r  rtahe  r 


«*  C 


(44) 
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u3u  .  v  3u  v  ^  1  3p 

rF*FJ?-i~*F  Sr  »  0  (45) 


u|^  +  I  |X  •>  Si  +  _i  |P  .  o 

Sr  r  «/6  r  pr  36 


(46) 


h 


+  (u 


2  .♦  V2, 


constant 


(4?) 


ui2.  +  v  £5.  ,  (48) 

3r  r  33  t 


Here  v.,  v  ars  the  velocity  components  in  the  r,  9  directions, 
respectively,  while 

e  =  0  for  plane  flow  (wedge  case) 

a  1  for  axisymnetr'ic  flow  (cone  case)  (49) 

If  the  shock  is  located  at  6  -  6  (r)  ,  the  shock  jun.p  conditions 
(10)  ~  i1!)  become 


VT  -s—  (r  sir,  3  i 
'■»  "dr  s 


V 


(r  ccs  8^'* 


5s  US  "  Vs* 


u  +  r  ? 1  v 


(50/ 

(51) 
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Ps  '  p-  laF*'  3ines»  „  i  , 

- - —  m  - - - 5-  (L  ~  —  ) 

P.V  1  +  *  1 


(52) 


Again  the  subscript  «•  denotes  proper-ties  upstream  of  the 
shock.  Note  that  the  shc:k  shape  6  (r)  is  to  be  determined  as 

5 

part  of  the  solution. 

Finally,  we  have  the  condition  of  no  flow  thru  the  body: 

v  =*  0  at  6  «  85  (53) 

where  0  *  9^,  a  constant,  locates  the  body. 

Were  dissociation  jbsent,  these  equations  would  have  the 
familiar  r-independent  iimilarity  solution.  However,  the 
equation  (48)  introducfs  a  characteristic  length  into  the  problem, 
and  the  solution  does  lepend  on  both  r  and  6  •  On  the  other 
hand,  the  problem  is  \yperbolic,  and  has  been  treated  success¬ 
fully  by  the  method  of  characteristics  (Capiaux  and  Washington 
1363;  see  also  Sedne.'  and  Gerber,  1963,.  1967). 

Two  versions  o)  the  series -expans ion  method  were  attempted 
for  these  problems.  The  first  is  what  Van  Dyke  (1966)  calls  the 
"method  of  series  ::*\mcation".  The  idea  here  is  to  expand  the 
variabiec  in  r,  eg., 

p  »  l  I!  *6)  r1  (54) 

i-P 

determin.-l^g  the  B-dependent  coefficients  like  p^(0)  by  solving 
(numer-Y-ally)  '.he  ordinary  differential  equations  derived  by 
subr-J-tutir.g  equations  like  (54)  into  the  governing  algebraic 
a.d  partial  differential  equations  and  then  collecting  terms  of 
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like  order  in  r.  Alternatively,  one  may  develop  the  variables 


in  double  series  in  r  and  8  ?  e.g., 

p  "  X)  pijrlr,j 

i,  j~o 

where,  bj  using 
n  = 


9-6,  (r) 


W 


(55) 


(56) 


as  an  independent  variable,  we  ease  the  satisfaction  of  the 
boundary  conditions  both  at  the  shocks  (n  *«  0)  and  at  the  body 


(n  =  1)  • 

In  either  approach,  the  hyperbolic  nature  of  the  problem 

is  reflected  in  the  determinability  of  the  coefficients  of  r1 

independent  of  those  of  x-5  for  j>i.  However,  as  in  equations 

„  i 

like  (22)-(24),  the  equations  governing  the  coefficients  of  r 
do  involve  the  coefficients  of  for  j<i,  so  that  the  coeffi¬ 
cients  must  be  determined  in  order,*  first  those  of  r^,  then 
those  of  r1,  etc.  In  the  full-expansion  method,  the 
coefficients  like  must  also  be  determined  in  order  of  their 
second  subscript  for  any  fixed  i ,  and  p ■ .  must  be  known  in 

X  '“i.  t  J 

order  to  find  p .  .  . 

In  both  methods,  the  coefficients  of  r1  are  found  by 


solving  a  boundary-value  problem,  in  which  the  shock  shape 
must  be  determined  so  as  to  produce  a  flow  field  which  meets 
the  body  boundary  condition  (53).  Specifically,  we  expand  the 


shock -shape  function 


(57) 


Then  0  ^  must  be  determined  by  iteration  (we  used  the  recant 
method)  so  that  the  coefficient  of  r1  in  the  expansion  of 
equation  ( h 3 )  vanishes. 

In  the  partial-expansion  approach,  this  boundary  condition 
becomes  simply 

yi(n)  “0  at  r,  -  1  (58) 

For  an  assumed  6^,  the  shock-jump  conditions  fix  the  values 
of  the  coefficients  like  and  at  n  13  0.  The  value  of 
vi(l)  for  that  was  then  determined  by  integrating,  from 
n  =  0  to  r,  1 ,  the.  system  of  ordinary  differentia)  equations 
formed  by  substituting  expansions  like  (54)  into  the  governing 
equations.  While  the  lower-order  coefficients  on  which  those  of 
r'1  depend  are  already  known  in  principle,  to  save  storage  they 
were  regenerated  by  solution  of  their  ordinary  differential 
equations  along  with  the  coefficients  cf  r1.  Of  course,  since 
8^  is  then  known  for  j<i,  the  problem  of  finding  the  lower -order 
coefficients  is  initial-value  in  nature.  Still,  the  partial- 
expansion  approach  was  found  to  be  much  more  time-consuming  than 
the  full -expansion  procedure,  and  no  more  accurate.  Perhaps 
the  time  problem  could  have  been  alleviated  by  storage  of  the 
known  lower-order  coefficients  (at  least  to  an  extent  which  would 
have  permitted  their  generation  by  interpolation) .  However, 
we  decided  to  spend  most  of  our  effort  or.  the  full-expansion 
method . 

The  derivation  of  the  recursion  formulas  for  the  series 
coefficients  in  the  full -expansion  method  proceeded  much  as  in 
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the  one-dimensional  problem  discussed  earlier.  To  illustrate, 
we  shall  consider  in  detail  the  treatment  of  the  differencial 
continuity  equation  (44),  wh..ch,  after  introducing  (56), 
becomes 

pii  +  up  +  (1+e)  (Qj^-Sp)  :u  -  epvz  +  v^j-  +  *  0  (60) 


Here  we  have  introduced  the  abbreviations 

*  =rt'Wlf +  : 

A  being  any  field  variable,  and 


(61) 


2  E  (es“0ij)  CGt  6 

=  (0g-eb)  cot  l(eb-e  jT,  +  eg]  (62) 


If  A  and  A  are  expanded  as-  in  equation  (55)  ,  on  substituting, 
such  expansions  into  equation  !6i)  ,  and  collecting  terms 
of  order  r1^ ,  we  find 


JL 

hi  ‘  X/  kek!i  Vm  '  li+l>  h-k.i+S 


k=0 


+  k  -  W 


(63) 


To  help  expand  (62),  we  r*ee:l  svroe  more  auxiliary  variables. 


Let 


a  =-  cos  (8s~0b)  d-n) 
s  2  sin  (0c,"Gi))  (1-n) 


and  expand  both  functions  in  series  ln'ke  (55).  Then 


L 

i=0 

v 1  "n 

lusi 


ciOr  *  008  <W 


iOr  "  sin  <W 


(  G  j  ) 


(65) 


Following  the  strategy  discussed  in  connection  with  equations 
(37)  and  (39) ,  we  differentiate  equations  (65)  with  respect 
to  r  and  collect  terras  of  order  r1  1  to  obtain 


ci0  “  cos(60-eb) 


if  i  -  0 


_T  1 


T  Ek<VVk(:>  "i- 


k=0  . 


si0  *  sin  lW 


.if  i  >  0 
if  i  »  0 


l  T  k  (9.-8.  6,  .)  c  ■  , 
x  L  k  b  kO  _-k 

k-0 


if  i  >  0 


(66) 


Similarly,  after  differentiating  (64)  with  respect  to  n  , 
we  find 


G5 


(6  V) 
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which  we  recast  as  a  continued  fraction  before  evaluation  at 
n  ~  1 .  In  the  case  of  flow  past  a  wedge,  this  recasting 
proved  to  be  unnecessary ,*  the  series  (69)  converged  quite 
rapidly  in  its.  raw  form.  The  cone  case  was  much  more 
difficult,  as  will  be  detailed  below. 

Thus,  in  essence,  the  full-expansion  method  simply  uses 
power-series  methods  to  integrate  the  ordinary  differential 
equations  which  govern  th»  9 -dependent  coefficients  of  the 
expansions  (54)  used  in  the  partial- expansion  method.  Gnce 
these  coefficients  are  available  at  a  fixed  value  of  9  (or 
rather  n  ),  a  rever  ting  of  the  resulting  series  in  r  into  a 
continued  fraction  in  r  was  necessary  to  obtain  final  results. 

Computations  have  been  carried  out  for  chi  four  cases 
treated  by  Capiat x  and  Washington  (196 i)  with  the  method  of 
characteristics,  see  Table  6.  Results  from  the  two  approaches 
are  compared  in  figures  2-4.  Where  they  disagree,  the  present 
results  are  presumed  to  be  correct,  at  least  for  relatively 
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small  values  of  x* ,  the  dimensionless  distance  along  the 
wedge  axis  {see  equation  (43) ) ,  our  justification  being  the 
convergence  of  our  series  expansions  evident  in  Table  7*.  The 
general  limitation  of  the  series  method  to  small  arguments, 
reflected  in  the  figure  by  plotting  the  difference  between 
11-  and  12-  term  results,  has  already  been  discussed.  The 
characteristics  method  is  not  inherently  so  limited,  but, 
being  an  initial-value  method,  it  may  become  too  costly  to 
pursue  the  calculations  to  completion..  This  is  especially 
true,  if,  as  in  the  Capiaux-Washington  calculations,,  one 
locates  the.  points  at  which  the  solution  is  to  be  determined 
at  the  intersections  of  characteristics  proceeding  downstream 
from  points  at  which  the  solution  is  known.  The  number  of 
points  at  which  calculations  are  made  then  increases  rapidly 
with  distance  downstream  of  the  wedge  apex.  It  would  have 
been  more  efficient  for  them  to  prescribe  the  solution 
points  a  priori,  using  Hartree’s  (1953)  interpolating  scheme 
to  integrate  the  compatibility  equations . 

?.s  it  is,  the  series  arid  characteristics  methods  are 
roughly  comparable  in  their  range  of  utility.  The  former 
method  is  much  moire  limi  ted  than  the  latter  in  case  4,  but 
does  do sc rise  the  entire  relaxation  rone  rather  well  even 
toon ;  it  >;  only  in  the  ct astant -property  equilibrium  region 
that  it  loses  accuracy.  in.  case  3 ,  the  present  method  is  clearly 

:  .  •  '.lota  '  hu .  in  th  i  -  case,  the  shock  wave  is  concave, 

.  .  ■  -  -  -’-r.  ‘  *\ ;  .  .  i  l '  J  .  i  .  ‘  .  ')  .  T  . .  t.  .  *  v  S 

r  os.il  .s  i:r  will  oh  they  converge  are  correct, 

.j  ..  r  -ivet  prog r by  substituting  the  series  results 

'■  *  .:!'■/  umvIoc  used. 
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equilibrium  value .  Capiaux  and  Washington  obtained  a  s  im.i  in  * 
result,  but  could  not  get  convergence  through  ref  in  orient  ol: 
their  mesh  size,  arid  so  discounted  the  over-shoot  as  a  numer¬ 
ical  inaccuracy  (to  be  sure,  the  overshoot  they  obtained 
was  much  more  severe  than  ours;  see  Figure  '/(c),  in  which 
Ax*  is  the  distance  between  the  first  two  points  on  the  wedge 
at  which  calculations  were  made,  and  so  is  indicative  of 
the  general  mesh  spacing) . 

Equations  (^4) -(52)  are  equally  applicable  to  flows 
past  wedges  and  cones;  the  only  difference  is  in  the  con¬ 
tinuity  equation  (44) ,  see  equation  (49) .  The  methods 
used  to  solve  those  equations  were  exactly  tne  same  for  both 
cases,  and  programs  were  written  in  which  the  input  value  of 
e  was  the  only  problem  identifier. 

Unfortunately ,  the  flow  Properties  a*.a  not  analytic  at 
the  cone  surface,  so  that  the  series  expansions  cannot  be 
made  to  converge  there.  To  see  this,  consider  equation 
(4S) ,  which,  under  the  transformation  (5b) ,  becomes 

ud  r  v~  -  CrJQ.-a,)?^  ( (p-pct)  e^V  1  -  c“a  Vbj )  (70) 

on  to  o 

The  fir"t  two  terms  in  the  partial  expansion  (in  r)  of  this 
equation  .ire 
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(73) 


Equation  (7  2)  shows  that,  as  would  be  expected,  the  flow  is 
"frozen”  (a  =  constant)  in  the  immediate  vicini'-v  of  the 
tip.  As  a  result,  x n  the  neroth-order  approximation  the 
flow  behaves  lik^  a  perfect  gas,  for  which  the  flow 
properties  are  certainly  analytic  in  9  and  n  . 

Thus,  using  (72)  ,  we  may  write  equation  (71)  as 

vo^  +  Wfo °i  -  £(">  <74) 


where  f  is  an  analytic  function  of  Tq,  p^,  and  aQ,  and  so 
is  an  analytic  function  of  as  well.  For  the  complementary 
solution  of  (74)  ,  we  try  a  Frobenius-type  expansion  of  the 
form 
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where*  the  t'orir.  of  uhe  oxponsien  for  v0  ref recta  the  body 
boundary  condition,,  Vq  «*  0  f  «  0.  Substituting  those 
expansion*  into  the  homogeneous  version  of  i/3),  wo  obtain 

0-  L  £  «ut(t+*)v0ij+1.i  +  (VV'o.j-h 

j-a  i«o  (77 

To  fix  6  ,  wo  insist  that  <*,-  •/  0 ,  and  so  obtain  the  ijo&icial 

X  w 

equation 


5  VCU  *  <VV  “00  - 0 


But,  collecting  terras  of  order  in  tha  continuity  equation 

(50) ,  we  get 

U+'H»b-Vu00  ♦  V01  “  0  (75) 


Jt  1 

5  *  TTT  <e 

Thus,  the  degree  of  dissociation  is  analytic  in  the  esse  of 
flow  past  a  wedge  (r*0),  but  haj  a  square -root  singularity 
at  the  surface  of  a  coic  . 

It  rs  rather  more  difficult  to  analyze  the  behavior  of 
the  other  flow  variables  at  the  body  surface.  If  we  aosuma 
expansions  like  that  j*  a,  in  (T<?  tor  &*.*  tSe  variables , 
we  find  p,  to  be  analytic  at  the  cone  surface,  ut  p  ^ , 
and  to  have  aguarj-roct  singularities  there.  Nc  result 
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for  vu  couin  be  proved.  liovever,  our  nur.er  J cal  calculations 
seerc  tc  indicate  -.hat  the  r -derivatives  of  p  and  v  are 
indeed  analytic  in  .  throughout  the  tlow,  but  that  the 
ether  variables  be  hi.ve  "irv-uiarr  tif  s  at  r~\  Our  evidence 
for  these  conclusions  is  exmoited  ^partially)  in  Ttfcle  £ , 
which  lists  the  coefficients  and  v^  ct'  our  original 
double-aeries  expansions  (55).  Since  the  v^'s  alternate  in. 
sign  and  decrease  in  magnitude  with  increasing  if  convergence 
of  the  series 

V  V,.",1 

i _  x  1 

1-0 

is  assured  at  least  for  g  <  1 .  However,  applying  Dorab's 
(1965)  procedures  to  tte  coefficients  a, . ,  we  obtain  the 
sequences  listed  in  Table  9  for  the  location  and  type  of  the 
singularity  closest  the  origin  of  c*.,  (n)  .  It  seams  clear  from 
these  data  that  this  singularity  is  at  n  =  i,  and  it  is  not 

toe  hard  tc-  believe  that  it  is  indeed  of  square -root  type. 

?  3 

The  coefficients  of  r  and  r'  in  tne  series  (54)  were 

alac  determined ,  and  behave  in  the  same  way  as  do  those  of 
1 

rA .  However,  the  higher  the  exponent  of  r ..  the  less  confidence 
we  have  in  our  determination  of  .  The  equations  for 
apparently  become  less  veil  conditioned  as  i  increases, 
raking  our  results  veil  infected  with  roundoff  errors,  so 
that  the  senes 


v.  1 

i  * 


while  not  diverging,  does  not  approach  a  fined  number  to 
any  acceptable  degree  of  accuracy. 

Sedney  <_nd  Gerber  (1567)  determined  the  first  two  terms 
of  an  expansion  of  the  flow  past  a  cone  of  a  gas  out  of 
vibrational  equilibrium,  and  s.iao  found  a  singular  behavior 
of  the  solution  at  the  cone  surface.  However,  their  expansions 
were  of  the  form 


P  -  P0U5  +  y  PX(C)  + 


where 


c  5  2*/y2 

*  is  a  stream  function. 
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—  *  ~CV  /  —  m  DV 
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y  y  sy 


x,  y  are  Cartesian  coordinates, 


(81) 


(82) 


(83) 


x  -  r  cos  e  ,  y  «  r  ain  8  (84) 

and  v  ,  v  are  the  velocity  components  in  the  x,  y  directions, 

y 

l/? 

It  was  not  clear  to  ua  .hat  the  c  singularity  of  their 

1/2 

solution  implied  a  (1-t»)  singularity  in  ours,  nor  is 
singular  behavior  evident  in  the  available  characteristics 
calculations  of  reacting  flow  past  cones  (Sedney  i  Gerber 
1963,  Spurck,  Gerber  and  Sedney  196’1)  .  It  would  be  of 
interest  to  attempt  the  removal  of  tiese  singularities, 
perhaps  by  matching  asymptotic  expann ions . 
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7 .  Steady  Supersonic  Flow  Past  a  Blunt  Body 

While  the  singularity  which  foiled  our  solution  of  the 
flow  past  a  cone  was  something  of  a  surprise,  the  existence 
of  a  singularity  ia  the  flew  field  of  the  corresponding 
blunt-body  problem  is  well  known  (Conti  and  Van  Dyke  1969, 
Vinokur  1970) .  (Specifically,  the  stagnation  point  of  such 
flows  is  singular.  This  phenomenon  is  connected  with  the 
fact  that,  since  the  residence  time  of  the  fluid  is  infinite 
at  the  stagnation  point,  the  flow  must  reach  equilibrium 
there,  no  matter  how  slow  is  the  reaction  rate. 

Nevertheless,  as  noted  in  the  Introduction,  one. of  the 
most  impressive  successes  of  the  series-expansion  method  in 
fluid  mechanics  has  been  its  solution  of  the  supersonic 
flow  of  a  perfect  gas  past  a  blunt  nedy.  Therefore  we 
wanted  to  see  how  it  would  perform  in  the  present  case  of  an 
ideal  dissociating  gas,  whatever,  a  priori  misgivings  we  may 
have  had.  Conti  and  Van  Dyke  (1969)  already  employed 
partial  expansions  of  the  form 

p  *  {r)  cos2  6  +  P2  (r)  sin20  +  O(sin^0) 

where  (r,9)  are  polar  coordinates  centered  at  the  center  of 
curvature  of  the  shock  ve»ve,  which  was  prescribed  as  either 
circular  (for  a  two-dimensional  calculation)  or  spherical 
(the  axisymmetric  case) .  They  did,  in  fact,  obtain  excellent 
results  for  the  flow  in  the  immediate  vicinity  of  the  stagnation 
point.  To  enable  approach  to  the.  singularity,  they  had  to 
introduce  as  an  independent  variable  (in  place  of  r)  the 


difference  on  the  axis  between  the  degree  of  dissociation 
and  its  equilibrium  value  (see  Conti  1966) . 

The  present  calculations  are  based  on  full  Taylor- 
series  expansions  of  the  form 


y* 

it 


Pi j (*"xs 


tr) ) S3 


(85) 


where  x»x  (z)  locates  the  bow  shock  wave  in  cylindrical 
coordinates  and  z  is  the  square  of  the  distance  from  the  axis 
of  symmetry.  Such  an  expansion  eases  the  satisfaction  of  the 
shock  jump  conditions  (10) -(14).  and  takes  advantage  of  the 
fact  that  p  is  even  in  the  cylindrical  radius  (to  insure 
that  all  independent  variables  behaved  similarly,  we  worked 
with  v//x  instead  of  the  radial  velocity  component  v,  which 
itself  is  odd  in  z.)  with 

y  i  x-xB(z) 

W  s  v/  /"z 

A  2  J*&  .  <*£  3A\ 
v5t  3y  / 

the  conservation  equations  (1)~(5)  become 


uf£  +  2wp  +  p|H  +  2pw  ♦  2ptf  »  0  (87) 

+  2wu  +  |  §£  °  0  (88) 

u|£  +  2ww  +  w2  +  L,  |2  «  o  (89) 

u|~  +  2vh  -  i  (u|£  +  2wp)  (90) 

u|y  +  2w5  *  3f-  (91) 
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The  shock  jump  conditions  (10)- (13)  may  bo  solved  for  the  flow 
properties  behind  the  bow  shock  as  follows: 
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( 93) 

(94) 

(95) 


Here 


Y  « 


sin6 
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(97) 


The  boundary  condition  of  no  flow  through  the  body  is  taken 
into  account  by  working  in  terms  of  a  stream  function  ii 
defined  by 


as  -  02  WZ  . 


ay 


a* 


dx 

0u  +  <!>  «  Oat  v*2**Q 

U.  i:  0.< 


(98) 


Then 

i  ,,  o  (99) 

z 

on  the  body . 

following  the  usual  approach ,  the  shock  shape 

x  \  i)  a;,  spec  it  ’  -?*d  'in  our  calculations ,  only  a  spherical 
«.*'  >•*.■!<  vavr  was  3 fuelled)  and  the  corresponding  body  shape 
do  urrTr.mv.;  an  parr,  of  the  solution  oy  searching  for  the 
v,o  ocs  of  vj/z.  The  recursion  relations  formed  by  substituting 
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the  series  like  (77)  into  the  governing  equations  turr  out 
to  be  solvable  in  sequence  of  the  total  order  i+ j  of  coeffi¬ 
cients  like  t>.  .  .  Amongst  coefficients  of  the  same  total 
''  -O 

order/  they  are  determined  in  order  of  the  first  subscript 
i.  Specifically,  the  snook  jump  conditions  (91) -  (94) 
determine  the  coefficients  with  subscript  i~0,  and  the 
differential  equations  (86) -(SO)  those  with  !>0. 

For  problems  in  more  than  one  independent  variable, 
there  is  no  procedure  directly  analogous  to  the  recasting  of 
a  one-variable  series  into  a  continued  fraction.  Thus,  to 
improve  the  convergence  of  two-dimensional  series  like  (84). 
we  must  somehow  put  them  into  one -dimensional  forms.  This 
can  be  done  in  any  number  of  ways.  Van  Tr.yl  (I960,  1967  , 

1971)  starts  v?ith  the  series  for  p  (e.g.1  in  the  form 

.  f  M-j 

p  3  L z3  Z  pij yl  (io° 

j»0  i=0 

recasts  the  polynomials  in  square  brackets  into  the  continued- 
fractions  P  (y)  (say),  and  then  calculates  p  from 


z*1  (y) 


(101) 


This  is  essentially  the  procedure  we  used  ir  our  analysis  of 
flows  past  wedges  and  cones,  as  described,  in  the  preceding 
section, except  that  the  series  corresponding  to  equation 
1100)  was  itself  recast  into  a  continued  fraction  in  what 
corresponded  to  the  present  variable  z . 
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However,  in  the  case  at  hand,  the  unavailability  of  the 
coefficients  p. .  of  total  order  i+j  beyond  some  finita  number 
(M  in  equation  (99))  means  that  the  coefficients  p..(y)  are 
known  with  decreasing  accuracy  as  j  increases .  Thus,  after 
experimentation  with  a  variety  of  procedures  similar  to 
Van  Tuyl’s  (Moran  1965),  we  found  it  preferable  to  rearrange 
series  like  (77)  into  the  form 


evaluate  the  polynomials  in  the  square  brackets  for  specified 
z/y,  and  then  recast  the  resultant  series  in  y  into  continued- 
fraction  form.  The  advantage  of  so  doing  is  that  the 
coefficients  of  y1  in  (101)  are  known  exactly  (within  roundoff 
errors)  up  to  the  point  at  which  the  series  in  y  is  truncated, 
since  the  coefficient  of  yx  consists  of  a  linear  combination 
of  the  coefficients  of  total  order  i.  Therefore,  our 
application  of  continued  fractions  may  be  interpreted  as 
looking  for  the  analytic  continuation  of  a  power  series  in  y 
on  a  curve  of  constant  z/y. 

Our  main  objective  was  to  find  the  shape  of  the  body  which 
generated  the  prescribed  (spherical)  shock  and  the  flow 
properties  on  that  body.  Once  its  series  coefficients  were 
available,  the  continued-fraction  expansion  of  fy/z  was  formed 
as  described  above  and  its  zeroes  determined  on  curves  of 
specified  z/y  by  the  secant  method. 

While  a  similar  procedure  worked  beautifully  in  the 
corresponding  perfect-oas  problem  (Moran  1970) ,  results  for 
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the  present  situation  wer*  very  poor.  Three  cases,  describe-;. 

in  Table  10,  were  examined  in  detail  for  comparison  with 

Conti  and  Van  Dyke's  (1969)  partial -expansion  results.  Ao 

shown  in  Figures  5,  in  the  two  cases  with  relatively  high 

rea  tion  rate,  wo  simply  were  unable  to  locate  the  body  with 

any  accuracy.  In  case  2,  the  present  rasuLts  give  no 

evidence  at  all  of  the  boundary  layer  at  the  stagnation 

point  in  which  the  fluid  makes  it*  final  approach  to  equilibrium. 

Only  in  case  3,  the  one  with  lowest  reaction  rate,  did  we 

get  any  indication  oi  the  presence  of  this  layer.  Also  in 

case  3,  our  results  for  the  degree  of  dissociation  seem  to 

suggest  the  presence  of  singularities  at  (and  beyond)  the 

stagnat  point,  While  this  is  correct  (Conti  and  Van  Dyke 

1969) ,  our  results  for  the  standoff  distance  (the  value  of 

y  at  which  y/z  •*  0  and  z  *  0)  appear  to  converge  to  an  answer 

different  from  Conti  ard  Van  Dyke's.  Wa  believe  their 

results  to  be  correct,  but  cannot  find  any  error  in  our 

calculations  either.  Our  program  was  checked  by  evaluating 

all  the  flow ‘properties  ind  their  derivatives  at  some 

point  within  the  shock  l*,yer  from  the  series  coefficients, 

and  determining  that  the  results  do  indeed  satisfy  the  governing 

equations. 

In  any  event,  the  apparent  failure  of  the  full  series- 
expansion  method  in  this  problem  was  not  entirely  unexpected. 

As  noted  above,  Conti  and  Van  Dyke  do  show  that  the  approach 
of  the  fluid  properties  to  their  etagnation  values  generally 
is  singular,  their  normal  derivatives  being  proportional  to 
some  negative  fractional  power  of  the  distance  from  the  body. 
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wjr  partial  success  in  describing  this  phenomenon  in  case  3 
is  just  as  surprising  as  the  apparent  convergence  in  that 
tase  to  the  wrong  result,  for  the  standoff  distance. 


8.  Diffraction  of  a  Plane  Shock  by  a  Blunt  Body 


Our  most  ambitious  project  in  this  series  was  the  deter¬ 
mination  of  the  time -dependent  flew  produced  when  a  plane 
ahock  wave  impinges  on  and  reflectt  from  a  solid  body  - 
specifically,  a  sphere  -  a  problem  cf  interest  in  connection 
with  the  blast-wave  loading  of  stationary  structures  and  with 
the  starting  of  shock  tub^a  and  tunnel Hero  we  break 
new  ground;  to  the  author's  knowledge,  this  problem  h»s  not 
been  solved  by  other  methods.  To  be  sure,  it  is  related  to 
the  well -studied  problem  o *:  finding  the  fl->w  about  a  body 
impulsively  accelerated  from  rest,  which  is  mainly  of  interest 
for  Its  final  steady  state. 

The  problem  is  illustrated  in  Figure  6.  Properties  are 

m 

known  in  the  undisturbed  region  1,  while  in  region  2,  behind 
the  incident  shock,  their  dependence  on  the  distance  from 
that  shock  can  be  found  by  the  procedure  described  *n  detail 
in  Section  2  of  this  report.  3oweyer  ,  for  all  caset 
considered  herein,  the  equilibrium  degree  of  dissociation  in 
region  2  was  less  than  10  6.  Since  the  other  flow  properties 
in  that  region  then  also  differ  negligibly  from  their 
equilibrium  values,  we  ignored  such  differences  in  our  analysis 
and  took  region  2  to  be  uniform. 

The  determination  of  the  flow  properties  in  region  3, 
between  the  body  and  the  reflected  ahock,  proceeds  much  as 
in  the  other  problems  described  in  this  repo. t,  with  the 
fallowing  major  differences » 
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as  on  two  space 


determined  simultaneously .  Specifically,  in  determining 
the  coefficients  of  total  or  dor  ifj+k.«n,  ail  iho«» 
coefficients  of  the  same  order  j  with  respect  to  the 
radial  coordinate  ^ust'  be  calculated  at  the  same  time, 
starting  with  j«b,  then  j**l,  etc.  The  equations  Involved 
are  linear  and  of  order  6(r»-j),  the  unknowns  being 
coefficients  at  the  series  expansions  of  p,  ur  pf  1/p, 
h,  end  T.  Coefficients  of  the  other  series  may  be 
determined  separately. 

The  shape  of  the  boundary  is  not  completely  known 

a  priori.  The  shock  shape  xr( z,t)  must  thus  b®  determined 

as  part  q£  th*i  solution,  partly  from  the  shock  jump 

conditions  and  partly  from  the  geometric  condition 

« 

that  the  incident  and  reflected  shocks  meet  the  body  at 

2 

the  sane  point,  whare  r  ■  Z^(t)  (say): 

vflt  m  *b**i*  *  *r  zl’  ^  *i05) 

Here  v#  is  the  spe-id  of  the  incident  shock. 

While  the  series  expansion  of  most  of  the  equations  is 
straight  forward,  ‘he  body  boundary  condition 

dx. 

v  *  w~3fJ’  at  x  -  '"'jjU)  (106) 

where  v  is  the  axi<;l  velocity  and  w  the  radial  valocify 
(divided  by  the  radius  to  produce  an  even  function  of 
the  radius) ,  la  relatively  troublesome .  The  problem  is  to 
evaluate  v  and  v  on  the  body  surface  in  terms  of  their 


series  coefficients.  Since,  icr  example. 

f 

vj  "1-  Vi^kvXfc(z)  “  {107; 

jx  -  xk-2)  ijk»0 


v**  nuat  introduce  the  auxiliary  functions 

YiU  t)  s  ixb(z)  •  2  ^ij**^* 

ijk=0 

Sines 

y.  -  (0(.t)  +  o/t))1 


%m  have 


1  ijk 


Thus 


if 


1+<<  i 


(103) 


(109) 

(110) 


j  k  1+  g 

v  -  E  *jtl<  E  E  E 

apoc^ix)  jk»0  £  “0  t.*0  i»G 

(111) 

The  easaa  for  which  confutations  were  carriod  out  are 
identified  in  Table  TI .  The  gas  properties  associated  with 
these  cases  are  about  the  ■same  as  in  \.apiaux  and  Washington’s 
(1963)  study.  The  value.',  chosen  for  the  ci.s  sedation  rate 
parameter  C  i.i  equation  *  7  /  correspond  roughly  to  body 
radii  of  1  era  {cases  1  a  a  3)  and  lb  or.  {case  2)  . 

Figures  7  and  3  show  the  distribution  of  pressure  and 
degree  of  disced aticn  along  the  body  surface  at  two  times 
early  in  the  process ,  while  chair  distributions  along  the 
axis  of  symmetry  sro  ahowr.  in  Figure*  9  and  10.  The 
convergence  i«  seen  to  be  slow  ev«a  when  v_t  ia  but  a  tenth 
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cf  the  sphere's  radius.  However,  rhe  otrertg  effect  of  .  .e 
dissociation  rate  parameter  C  is  quite  clear.  In  th  • 
corresponding  perfect-gas  problem,  the  raethod  performed  much 
better,  yielding  useful  results  in  the  vicinity  of  the 
stagnation  point  for  vst/r^  up  to  2.0  or  i^ore  iHoran  and 
Van  ifoorheir.  19d9)  .  Of  course,  it  wes  to  be  expected  that 
the  present  problea,  with  its  relaxation  process,  would  be 
such  ;aore  difficult  than  the  other. 
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Figure  Captions 


Figure  1.  Temperature  behind  normal  shock  in  Llghthill  gas 
(normalized  by  \£/R) .  Cases  identified  in  Table  2.  Solid  lines 
arc  Freeman's  (1958)  results;  circles,  continued-fraction  results 
{34  and  33  terms)  . 

Figure  2.  Angle  between  tangent  to  shock  attached  to  wedge  in 
Lighthill  gaa  and  freestream  direction  (radius) .  Solid  lines 
are  Capiaux  and  Washington1 &  (1963)  results;  circles,  continued- 
fraction  results  (11  and  12  terms),  (a)  Case  1  of  Table  6. 

Figure  2 (b) -  Casa  2  of  Table  6 . 

Figure  2(c).  Case  3  of  Table  6. 

Figure  3.  Presume  distribution  on  wedge  in  Llghthill  gas 

2 

(normalized  by  pjrJ').  For  legend  see  Figure  2.  (a).  Case  1 

of  Tabic  6. 

Figure  i (b\  Case  2  of  Table  6. 

Figure  2(c)  Oise  1  of  Table  6. 

Fignre  4,  Temperature  distribution  cr,  vodge  in  Lighthill  gaa 
(normalized  by  v^/R)  .  For  legend  see  Figure  '2,  (a).  Case  1 

of  Table  6 . 

g-'-r,.  5  (b>  .  Case  or  Tafo3e  6. 

-  ■  i  .  e.  4  (c)  ,  ;  ase  4  of  Tide  6. 


G2 


Figure  5.  Streamfunction  divided  by  Bquare  oil  radiur  (i{//z)  and 
degree  of  dissociation  (a)  on  axis  of  symmetry  in  shock 
layer  of  blunt  body  which  supports  spherical  shock  wave  in 
Lighthiii  pas.  Y  its  distance  downstream  of  nose  of 
shock  (normalized  by  shock  radius) . 

Solid  lines  are  Conti  and  Van  Dyke's  (1969)  results;  circles, 
continued-fraction  results  (34  and  3 5  terms) . 

(a).  Case  1  of  Table  10. 

Figure  5(b).  Case  2  of  Table  10. 

Figure  5  (c) .  Case  3  of  Table  10  > 

Figure  6.  Diffraction  of  plana  shock  by  solid  body. 

Figure  7.  Pressure  distribution  on  surface  of  sphere  of  radius  1 

2 

after  impingement  of  shook  (normalised  by  p-Vj  ) .  Times 
normalised  by  l/vy .  Continued*  fraction  results  (10  and  11  terma) . 
Circles,  case  1  of  Table  XX;  triangles,  case  2;  squares,  case  3. 

Figure  8.  Degree  of  dissociation  on  surface  of  sphere  of  radius  l 
after  impingement  of  shock.  For  legend  see  Figure  7. 

Figure  9*  Pressure  distribution  of  axis  between  reflected  shook 

2 

end  sphere  of  radius  (normalised  by  PjV3  ) .  For  legend  see 
Figure  7. 

Figure  10.  Degree  of  dissociation  on  axis  between  reflected  shock 
and  sphere  of  radius  1.  For  legend  see  Figure  7. 
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